Oh 

O 
43 



DO-TH 95/13 
RAL-TR-95-042 
August 1995 
updated version (December 1995) 



Next-to-Leading Order Radiative Parton Model Analysis 
of Polarized Deep Inelastic Lepton Nucleon Scattering 

> 

M. Gliick, E. Reya, M. Stratmann 

oo . 

Universitiit Dortmund, Inst it ut fiir Physik, 

in ■ 

Q\ ' D-44221 Dortmund, Germany 



P-i' 

D ! W. Vogelsang 



Rutherford Appleton Laboratory 
X ' Chilton, Didcot, Oxon 0X11 OQX, England 



Abstract 



A next-to-leading order QCD analysis of spin asymmetries and structure functions in 
polarized deep inelastic lepton nucleon scattering is presented within the framework of 
the radiative parton model. A consistent NLO formulation of the Q 2 -evolution of polarized 
parton distributions yields two sets of plausible NLO spin dependent parton distributions 
in the conventional MS factorization scheme. They respect the fundamental positivity 
constraints down to the low resolution scale Q 2 = ^ NLO = 0.34 GeV 2 . The Q 2 -dependence 
of the spin asymmetries A?' n ' (x, Q 2 ) is similar to the leading-order (LO) one in the range 
1 < Q 2 < 20 GeV 2 and is shown to be non- negligible for x-values relevant for the analysis 
of the present data and possibly forthcoming data at HERA. 



1 Introduction 



Recently, a leading order (LO) QCD analysis of polarized deep inelastic lepton nucleon 
scattering has been performed IJ within the framework of the radiative parton model. 
The first moments Af(Q 2 ) of polarized parton distributions Sf(x,Q 2 ), 

Af(Q 2 ) = f dx5f(x,Q 2 ) , (1.1) 
Jo 

where f — u,u, d, d, s, s, g, were subject to two very different sets of theoretical constraints 
related to two different views concerning the flavor SU (3) [577 (3) /] symmetry properties of 
hyperon /3-decays. One set ('standard' scenario) assumed an unbroken 577(3)/ symmetry 
between the relevant matrix elements while the other set ('valence' scenario) assumed an 
extremely broken 577(3) / symmetry reflected in relating M the hyperon /3-decay matrix 
elements to the first moments of the corresponding valence distributions. Both scenarios 
gave satisfactory descriptions [|IJ of the measured spin-asymmetries [3-9] A\' n {x ) Q 2 ) ~ 
gi ,n (x, Q 2 )/ Fi' n (x, Q 2 ), although the polarized gluon density 5g(x,Q 2 ), which enters in 
LO only via the (Revolution equations, was only weakly constrained by present data. 
The total helicity carried by quarks 

A£(Q 2 ) = £ {Aq{Q 2 ) + Aq{Q 2 )) , (1.2) 

q=u,d,s 

which is Q 2 -independent in LO, turned out to be AS ~ 0.3 in both scenarios with an 
average total gluonic helicity Ag(Q 2 = 4GeV 2 ) ~ 1.5. A specific feature of our radiative 
LO analysis is that the polarized leading twist parton densities Sf(x, Q 2 ) are valid down 
to Q 2 = fj 2 ^ ~ 0.23 GeV 2 and that the fundamental positivity constraints 

Sf(x,Q 2 )\ </(x,Q 2 ) (1.3) 

are respected down to this low resolution scale Q 2 = /i| as well. The 'standard' scenario 
requires a finite total strange sea helicity of As = As ~ —0.05 in order to account for a 
reduction of T{ with respect to the Gourdin and Ellis and Jaffe estimate [1(| 

i> + ^ + 3lr 



n,Ej = T^( F + D ) + ^(3F -D)~ 0.185 (1.4) 
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where 

T\{Q 2 ) = j l dxgl{x,Q 2 ) . (1.5) 
J o 

Within the 'valence' scenario, on the contrary, a negative light sea helicity Am = Ad 
At I ~ —0.07 suffices (As = As = 0) for reducing T\ EJ . Furthermore, in both scenarios 
we Q predict T\ ~ 0.15 and ~ —0.06 in agreement with recent experiments [5-9], with 
the Bj0rken sum rule being manifestly satisfied. 

While our LO analysis was being completed, a full next-to-leading order (NLO) calcu- 
lation of all polarized two-loop splitting functions 8P[p(x), i,j = q, g, in the conventional 
MS factorization scheme has appeared for the first time ||11|| . It is the purpose of this arti- 
cle to present first a consistent NLO formulation of spin-dependent parton distributions, 
making use of the NLO results of ref.[|TTJ, in particular for (Mellin) n-moments of struc- 



ture functions and parton densities where the solutions of the NLO evolution equations 
can be obtained analytically. Using these formal results we then proceed to perform a 
quantitative NLO analysis of Ai' n (x, Q 2 ) and gi ,n (x, Q 2 ), and will present two sets of NLO 
Sf(x, Q 2 ) for the two scenarios discussed at the beginning. Since most NLO analyses con- 
cerning unpolarized hard processes and parton distributions have been performed in the 
MS factorization scheme, it is convenient to remain within this factorization scheme also 
for polarized hard processes and spin- dependent parton distributions. This is particularly 
relevant for the parton distributions which have to satisfy the fundamental positivity 
constraints (1.3) at any value of x and scale Q 2 , as calculated by the unpolarized and 
polarized evolution equations, within the same factorization scheme. In addition we also 
repeat our previous LO analysis [1J since new data have been published very recently 

EE- 

2 NLO Parton Distributions and their Q 2 -Evolution 

Measurements of polarized deep inelastic lepton nucleon scattering yield direct informa- 
tion [3-9,12] on the spin-asymmetry 

a n ( n 2\ _ 9ii x ,Q 2 ) _ 9i{x,Q 2 ) (r> 1 s 

1 { ' V } - F?(x, Q') F»(x, Q*)/ [2x(l + W(x, Q 2 ))) ' [ZA) 
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N = p,n and d = (p+n)/2, where in the latter case we have used gf = \{g\+gi )(1 — |^d) 
with ojb = 0.058 0, |!|; R = Fpj2xF\ = (F 2 — 2xF\) /2xF\ and subdominant contributions 
have, as usual, been neglected. In NLO, A±(x, Q 2 ) is related to the polarized (Sf N ) and 
unpolarized (f N ) quark and gluon distributions in the following way: 



g?{x,Q 2 



{Sq N (x,Q 2 ) + Sq N (x,Q 2 ) + 



+ 



a s (Q 2 ) 



2tt 



5C q * (5q N + 5q N ) + j5C g * Sg 



(2.2) 



with the convolutions being defined by 



(C*q)(x,Q 2 )= f d JL C (*) q(y,Q 2 ) 

J* v \yj 



(2.3) 



and where the appropriate spin- dependent Wilson coefficients in the MS scheme are given 



by (see [11], for example, and references therein) 

2 /ln(l — x) \ 3 



SC q (x) 



l + x< 



+ 2 + x 



SC g (x) 



2T f 



(2x 



\ 1 — x 



1 i ( In ^ ^ — 1 

x 



2(1 -x) 



1 + x 2 
1 — X 



In x + 



2(1 -x) 



(2.4) 
(2.5) 



with Cp = 4/3 and Tf = f/2. Here / denotes, as usual, the number of active flavors 
(/ = 3). The NLO expression for the unpolarized (spin-averaged) structure function 

F N 

(x,Q 2 ) is similar to the one in (2.2) with 5f(x,Q 2 ) — > f(x,Q 2 ) and the unpolarized 
Wilson coefficients are given, for example, in |13| . Henceforth we shall, as always, use the 
notation 5q p = 5q and q p = q. Furthermore the NLO running coupling is given by 



a s (Q 2 



(3, lnlnQ 2 /A 



2 

MS 



4vr /? lng 2 /A|s Pi (lnQVA 



2 

MS 



(2.6) 



with /3b = 11 - 2//3, A = 102 - 38//3 and being given by [|14 



A. 



(M,5) 

MS 



248, 200, 131 MeV . 



The number of active flavors / in a s (Q 2 ) was fixed by the number of quarks with m 2 < 
Q 2 taking m c = 1.5 GeV and wif, = 4.5 GeV. The marginal charm contribution to g± , 
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stemming from the subprocess 7*g — > cc | L3(| , will be disregarded throughout. The charm 
contribution to is also small in the kinematic range covered by present polarization 
experiments. 

For calculating the NLO evolutions of the spin-dependent parton distributions Sf(x, Q 2 ) 
in (2.2) we have used the well known analytic NLO solutions in Mellin n-moment space 

with the n-th moment being defined by 
»i 



(see, e.g., refs.[13, 16 



5f n (Q 2 ) = / dxx n ~ l 5f{x,Q 2 ) . 
Jo 

These (Revolutions are governed by the anomalous dimensions^ 

17 = ±1 
i,j =q,9 



(2.7) 



hi 



a 



Air 



5^ )n + ( 



Air 
Air 



(2.8) 
(2.9) 



whose detailed n-dependence will be specified in the Appendix. The non-singlet (NS) 
parton densities evolve according to [|13|, [RJ 



Wns v (Q 2 ) 



Air 



2/3 



' a s (Q 2 ) 

MQl) 



INS r](Ql) 

(2.10) 



with the input scale Qq = //jvlo = 0.34 GeV 2 referring to the radiative \TA~\ NLO input 
(Unlo = Hho) to be discussed later. Furthermore, opposite to the situation of unpolarized 
(spin-averaged) parton distributions [IB], 5qNSr]=+i corresponds to the NS combinations 
5u — 5u = buy and 5d — 5d = 5dy, while SqNSri=-i corresponds to the combinations 
5q + 5q appearing in the NS expressions 



Sq 3 = (5u + Su) - (5d + 5d) , Sq 8 = (5u + 8u) + [Sd + 5d) -2(5s + 5s) 



(2.11) 



It should be noted that the first (n = 1) moments 8q NS - = Aq NS _ of these latter 
SU(3)f diagonal flavor non-singlet combinations do not renormalize, i.e. are indepen- 
dent of Q 2 , due to the conservation of the flavor non-singlet axial vector current, i.e. 



1 Alternatively one can of course use instead the LO and NLO splitting functions SP^ n = —Sj^^/A 



and 5R 



(l)n 



-(?7-^ n /8, respectively (see, e.g., ref.Jl7t). 
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&lqq l = $1ns(V = ~ 1) = ( see Appendix). The evolution in the flavor singlet sector, 



iqq 

i.e. of 



5T, n {Q 2 ) = Y, [bq n {Q 2 ) + b<f{Q< 



(2.12) 



q=u,i 



and Sg n (Q 2 ) is governed by the anomalous dimension 2x2 matrix (2.9) with the explicit 



solution given by eq.(2.9) of ref . [ 16 1 where 7 — > £7 as given in the Appendix. 



Having obtained the analytic NLO solutions for the moments of parton densities, 
Sf n (Q 2 ), it is simple to (numerically) Mellin-invert them to Bj0rken-a; space as described, 
for example, in |16| or ||17| . The so obtained 5f(x,Q 2 ) have to be convoluted with the 
Wilson coefficients in (2.2) to yield the desired gi(x,Q 2 ). Alternatively, one could insert 
Sf n (Q 2 ) directly into the n-th moment of eq.(2.2), 



+ 



2tt 



with the moments of (2.4) and (2.5) given by 



SC™ = C F 



S 2 (n) + (S 1 (n)) 2 + ri 



1 



2 n(n+l, 
" _J {S 1 {n) + l)-± + 



Si(n)+ -r + — + 



f2.2'1 



n 2 2n n + 1 



'2.4' 



(2.5') 



n(n+l)' " v ' ' n 2 n(n + l) 

with Sk(n) defined in the Appendix. The full expression (2.2') can now be directly (nu- 
merically) Mellin-inverted [JT6], [XT] to yield gi(x, Q 2 ) without having to calculate any con- 
volution (2.3). 

The LO results are of course entailed in all these expressions given above, by simply 
dropping all the obvious higher order terms (/?i, S^ 1 ', SC q>g ) in all relevant equations 
stated above. 



It should be noted that the first (n = 1) moment Ti(Q 2 ) = g\(Q 2 ) in (1.5) is, according 
to (2.2'), simply given by 



Ti(Q 2 



« S (Q 2 



7T 



Aq(Q 2 ) + Aq(Q 2 



(2.13) 
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where we have used the definition (1.1) and 8Cg = — 3CV/2 and 8Cg = according to 
(2.4') and (2.5'), respectively. Thus, the total gluon helicity Ag(Q 2 ) does not directly 
couple to ri(Q 2 ) due to the vanishing of the integrated gluonic coefficient function in the 
MS factorization scheme. This vanishing of AC g = SCg, which has been some matter of 
dispute during the past years (for reviews see, for example, [[L2], [18], |19[) originates from 



the last term in (2.5) proportional to 2(1 — x). Since this term derives from the soft 
non-perturbative collinear region |2(J where k\ ~ m 2 <C A 2 , it has been suggested [18, 
21-23] to absorb it into the definition of the light (non-perturbative) input (anti) quark 
distributions 5 Q (x,Q 2 = Qq). This implies that, instead of 5C g (x) in (2.5), one has 

5C g (x) = 2T f {2x- 1) - l) (2.14) 



which refers to some different factorization scheme |18, [22], ^] with the n-th moment 
given by 

SC n n = 5C n n - 2T f —^ — - (2.14') 
9 9 } n(n + l) v ; 

and 5Cg given in (2.5'). Thus AC g = 5Cg = — 2Tf and Ag(Q 2 ) would couple directly 

21 1 to Txi^Q 2 ) in (2.13) via — (a s /67r)Ag(Q 2 ) according to the gluonic term in the curly 



brackets of (2.2') for / = 3 flavors. Therefore the gluonic contribution on its own could 
account for a reduction fL8 , 21, ^] of the Ellis- Jaffe estimate (1.4), as required by exper- 



iment, without the need of a sizeable negative total (strange) sea helicity as discussed in 
the Introduction. One could of course choose to work within this particular factorization 
scheme or any other scheme. In this case, however, one has for consistency reasons to 
calculate all polarized NLO quantities (8C*,8j!jp n , etc.), and not just their first (n = 1) 
moments, in these specific schemes as well as also NLO subprocesses of purely hadronic 
reactions to which the NLO parton distributions are applied to. The transformation^ 



SCg — > 8Cg in (2.14') implies of course also a corresponding modification [13, [24j] of the 
NLO anomalous dimensions 87^™ — > 8^\^ n . 



2 This transformation cannot even be uniquely defined for the spin dependent case. This is partly in 
contrast to the unpolarized situation [|l3| [24], [ll| where the energy-momentum conservation constraint 
(for n = 2) is used together with the assumption of its analyticity in n. 
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3 Quantitative NLO Analysis 



In fixing the polarized NLO input parton distributions Sf(x, Q 2 = /-tjVLo) we f° now closely 
our recent LO analysis |]]. We still prefer to work with the directly measured asymmetry 
Ai(x, Q 2 ) in (2.1), rather than with the derived gx(x, Q 2 ), since possible non-perturbative 
(higher twist) contributions are expected to partly cancel in the ratio of structure functions 
appearing in A±(x, Q 2 ), in contrast to the situation for g^(x,Q 2 ). Therefore we shall 
use all presently available data [4-9] in the small-x region where Q 2 ^ 1 GeV 2 without 
bothering about lower cuts in Q 2 usually introduced in order to avoid possible higher twist 
effects as mandatory for analyzing gi(x, Q 2 ) in the low-Q 2 region. The analysis affords 
some well established set of unpolarized NLO parton distributions f(x, Q 2 ) for calculating 
(x, Q 2 ) in (2.1) which will be adopted from ref. |14j, i.e. our recent updated NLO (MS) 



dynamical distributions valid down to the radiative input scale Q 2 = [i 2 NLO — 

0.34 GeV . 

The searched for polarized NLO (as well as LO) parton distributions 5f(x,Q 2 ), com- 
patible with present data [4-9] on A±(x, Q 2 ), are constrained by the positivity require- 
ments (1.3) and for the SU(3)f symmetric 'standard' scenario by the sum rules 

Au + Au - Ad - Ad = g A = F + D = 1.2573 ± 0.0028 (3.1) 
Au + Aw + Ad + Ad - 2(As + As) = 3F - D = 0.579 ± 0.025 (3.2) 

with the first moment Af defined in (1.1) and the values of qa and 3F — D taken from 
Z5 |. It should be remembered that the first moments Ag 3)8 of the flavor non-singlet 



combinations (2.11) which appear in (3.1) and (3.2) are Q 2 -independent also in NLO due 
to S^ns (rj = — 1) = 0, according to eq.(A.9). 

As a plausible alternative to the full SU(3)f symmetry between charged weak and 
neutral axial currents required for deriving the 'standard' constraints (3.1) and (3.2), we 
consider a 'valence' scenario |T], ^] where this flavor symmetry is broken and which is 
based on the assumption [0 that the flavor changing hyperon (3-decay data fix only the 
total helicity of valence quarks Aq v (Q 2 ) = Aq — Aq : 

AM/4ao) - ^M^nlo) = 9a = F + D = 1.2573 ± 0.0028 (3.1') 



AM/4lo) + &Mt*NLo) = 3F — D = 0.579 ± 0.025 . (3.2') 

Although at the input scale Au(n 2 NLO ) = Ad(n% LO ), isospin symmetry will be (marginally) 
broken by the NLO evolution, i.e. Au(Q 2 ) ^ Ad(Q 2 ) for Q 2 > [i% L0 . In addition 
we shall assume a maximally SU(3)f broken polarized strange sea input 5s(x, [i 2 NLO ) = 
Ss(x, njyLo) = in our 'valence' scenario, which in addition is compatible with the SU(3) / 
broken unpolarized radiative input s{x ) ^i 2 NLO ) = of ref.[14]]. Such a choice is feasible in 
the 'valence' scenario since, due to eq.(2.13) and (2.11), we have in general 



a s (Q 2 



7T 



(3.3) 



Since Auy(Q 2 ) — Ady(Q 2 ) and Auy(Q 2 ) + Ady(Q 2 ) decrease only marginally with Q 2 as 
compared to their input values in (3.1') and (3.2') [cf. Table 1' below], and furthermore the 
dynamical isospin breaking Am ^ Ad at Q 2 > h%lq is small, eq.(3.3) can be approximated 
by 



r?' n (Q 2 



± —IF + D) + —(3F -D) + 
12 v ' 36 V ' 



+ — (10Ag(Q 2 ) + As(Q 2 ) + As(Q 2 



a s {Q* 



71 



(3.4) 



Therefore a light polarized sea Aq = (Au + Ad)/ 2 < can account for the reduction of 
the Ellis- Jaffe estimate (1.4) for Y\{Q 2 ), say, as required by recent experiments [|], |J. 



In the 'standard' scenario we need, on the contrary, a finite sizeable As(Q 2 ) < since, 
due to eq.(2.13) and (2.11), 

r?' n (Q 2 ) = [±^Ag 3 + l 6 Aq 8 + \ (As(Q 2 ) + As(Q 2 ))] (l - (3.5) 

with the Q 2 -independent flavor non-singlet combinations Aq^s being entirely fixed by 
eqs.(3.1) and (3.2); for the singlet combination in (1.2) we used A£(Q 2 ) = Ag 8 + 
3 (As(Q 2 ) + As(Q 2 )). It should be noted that, in contrast to the LO [JIJ finite 
As(Q 2 ) will be generated dynamically in NLO for Q 2 > n% LO even for a vanishing input 
As(fi] VLO ) = due to the non- vanishing NLO A7W in (A. 9): The resulting 'dynamical' 
As(Q 2 ) < is about an order of magnitude too small to comply with recent experiments 
|5|, H which typically yield, for example, T^Q 2 = 3 GeV 2 ) ~ 0.12 — 0.13, i.e. sizeably 



smaller than the naive estimate (1.4). We therefore have to implement a finite strange 
sea input As(fi 2 NLO ) = As(fi 2 NLO ) < for the 'standard' scenario, in order to arrive at 
As(Q 2 = 3 — 10 GeV 2 ) ~ —0.05 as required [IJ by recent experiments. 

Apart from applying the above scenarios for the polarized input distributions to 
A^(x,Q 2 ) rather than to gi(x,Q 2 ), the main ingredient of our NLO analysis is the 
implementation of the positivity constraints (1.3) down to [14] Q 2 = fi% LO — 

0.34 GeV 



(and to Q 2 = n 2 LO = 0.23 GeV 2 in LO) which is not guaranteed in the usual (LO) studies 
done so far (recently, e.g., in [26-29]) restricted to Q 2 > Ql = 1 — 4 GeV 2 . We follow here 
the radiative (dynamical) concept which resulted in the successful small- x predictions of 
unpolarized parton distributions as measured at HERA |TJ|, [H], |3(| . A further advantage 



of this analysis is the possibility to study the Q 2 -dependence of A±(x, Q 2 ) in the small- x 
region over a wide range of Q 2 || which might be also relevant for the forthcoming polar- 
ized experiments at HERA. In addition it will be important to learn about the reliability 
of perturbative calculations by comparing the LO with the NLO results; a reasonable 
perturbative stability of all radiative model predictions will be indeed observed for mea- 
surable quantities such as Af(x, Q 2 ) and (x, Q 2 ), as is the case for spin-averaged deep 
inelastic structure functions |T^, |3T)|j . 

Turning to the determination of the polarized NLO (LO) parton distributions Sf(x, Q 2 ) 
it is helpful to consider some reasonable theoretical constraints concerning the sea and 
gluon distributions, in particular in the relevant small- x region where only rather scarce 
data exist at present, such as color coherence of gluon couplings at x ~ and helicity 
retention properties of valence densities as x — > 1 | 3~If . We follow here very closely the 



procedure and ansatze of ref.|IJ. Subject to these constraints we employ the following 
general ansatz for the NLO (LO) polarized parton distributions: 

5q v (x,fi 2 ) = N qv x a «v q v (x,fi 2 ) 
5q(x, fx 2 ) = NqX a " (1 - x) b " q(x, /i 2 ) 
5s(x, fi 2 ) = 5s(x, fi 2 ) = N s Sq(x, /i 2 ) 

5g(x,fx 2 ) = N g x a9 (1 -x) bB g(x,fx 2 ) (3.6) 
with the NLO (LO) unpolarized input densities being taken from ref.[H] and, for obvious 



reasons, we have not taken into account any SU(2)f breaking input (5u 7^ 5d) as is 
apparent from our ansatz for 5q = 8u = Sd proportional to q = (u + d)/2 which should be 
considered as the reference light sea distribution for the positivity requirement (1.3). In 
the 'standard' scenario our optimal NLO densities at Q 2 = /U^ i0 = 0.34 GeV 2 are given 
by 

N Uv = 0.6586 , a Uv = 0.18 ; N dv = -0.3392 , a dv = 
Ng = -0.525 , a q = 0.51 , h- q = ; N s = 1 

N g = 10.47 , a g = l.l , 6 9 = 4.3 (3.7) 

corresponding to y 2 = 104.5/125 d.o.f. and respecting eqs.(3.1) and (3.2) which are the 
basis of almost all analyses performed so far. 3,4 Since new polarized p and d data have 
appeared very recently H |7|, we have also again determined the LO distributions in the 
'standard' scenario at Q 2 = jjL 2 LO = 0.23 GeV 2 : 

N Uv = 0.6392 , a Uv = 0.16 ; N dv = -0.3392 , a dv = 
Ng = -1.348 , a q = 0.93 , b- q = 0.09 ; N s = 1 

N g = 11.8 , a g = 0.76 , b g = 6.84 (3.7') 

corresponding to x 2 = 108.0/125 d.o.f. and which supercedes our LO input fit of ref . ]l[] . 
The fact that Ss(x, jj 2 ) 7^ in (3.7) and (3.7'), i.e. N s = 1, contradicts somewhat our 
purely radiative input |L4j s(x,/i 2 ) = s(x,/i 2 ) = 0, but for Q 2 ^ 1 GeV 2 the positivity 



inequality (1.3) is already satisfied, in particular at large values of x. In this respect the 
input for the 'valence' scenario with the extreme 577(3)/ breaking ansatz Ss(x, fj, 2 ) = 
is more consistent as far as our radiative (dynamical) approach is concerned: For the 
577(3)/ broken 'valence' scenario, based on the constraints (3.1') and (3.2'), our optimal 

3 It is interesting to note that, within our radiative approach with its longer Q 2 -evolution 'distance', 
a finite (negative) strange sea input 5s(x, /ijvLo) ^ s always required by present data even if one uses 5C g 
in (2.14) or (2.14') [21-23]. For the latter case, 8s has to be at least half as large as in (3.5) which is 
based on the MS SC g in (2.5) or (2.5'). This holds true even for a maximally saturated input gluon 
[Sg(x, nff LO ) = g(x, /Ujvlo)] ^° ^ e discussed below. 

4 It should be noted that our fit result N s = 1 in (3.7) implies via (3.6) an 5/7(3)/ symmetric sea 
input. In NLO at Q 2 > h% L q this symmetry becomes dynamically broken via the Q 2 - evolution. In view 
of present scarce data such SU (3) / [and SU (2) A breaking effects are entirely negligible for quantitative 
analyses. 
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NLO input corresponds to the following parameters in (3.6) at ^ 2 NLO — 

0.34 GeV 2 : 

N Uv = 0.6708 , a Uv =0.19 ; N dy = -0.3693 , a dv = 0.03 
Nq = -0.6 , a q = 0.5 , bg = ; N s = 

N g = 9.87 , a g = 1.05 , b g = AAA (3.8) 

corresponding to x 2 — 104.6/125 d.o.f. Similarly our new LO input at fi 2 LO = 0.23 GeV 2 
in the 'valence' scenario is given by 

N Uv = 0.6272 , a uv = 0.15 ; N dy = -0.3392 , a dv = 
N q = -1.455 , ag = 0.88 , b- q = 0.13 ; N s = 

N g = 7.4 , a g = 0.6 , 6 9 = 5.93 (3.8') 

corresponding to \ 2 = 108.6/125 d.o.f. and which supercedes our LO input fit of ref.Q. 
Finally, similarly agreeable NLO 'valence' scenario fits to all present asymmetry data 
shown below (with a total x 2 of 106.2 to 107.6 for 125 data points) can be also obtained for 
a fully saturated (inequality (1.3)) gluon input Sg(x, f^ 2 NLO ) = g{ x i Pnlo) as we ^ as ^ or ^ ne 
less saturated Sg(x, [i 2 NLO ) — %g{%, I^nlo)- ^ purely dynamical |32[ input 8g(x, H 2 NLO ) — 
is also compatible with present data, but such a choice seems to be unlikely in view of 
Sq(x, fi 2 NLO ) 7^ 0; it furthermore results in an unphysically steep |32| 5g(x,Q 2 > H 2 NLO ), 



being mainly concentrated in the very small-x region x < 0.01, as in the corresponding 
case [|L6"|, j33fl for the unpolarized parton distributions in disagreement with experiment. 
Similar remarks hold also for a LO analysis as well as for the 'standard' scenario. 

A comparison of our results with the data on A±(x, Q 2 ) is presented in Fig.l. The 
LO and NLO results in the 'standard' scenario are perturbatively stable and almost 
indistinguishable. The same holds true for the results in the 'valence' scenario which are 
not shown separately since they almost coincide with the 'standard' ones in Fig.l. As 
already mentioned, fit results using a 'saturated' gluon 5 g = g or Sg = xg, or even 5g = 
at Q 2 = fJ-io nlo are ver y similar to the ones shown in Fig.l. Note that A±(x,Q 2 ) — > 
const, as a: — > 1. The recently measured Q 2 - dependence of A^(x, Q 2 ) is compared with 
our theoretical results down to Q 2 = 0.4 GeV 2 in Fig.2. The difference between our LO and 
NLO results in the small-Q 2 region is mainly due to different LO [p? LO = 0.23 GeV 2 ) and 
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NLO (fj? NLO = 0.34 GeV 2 ) input scales. The more detailed Q 2 -dependence of A^(x, Q 2 ) is 
presented in Fig. 3 for some typical fixed x values for 1 ^ Q 2 < 20GeV 2 relevant for present 
experiments. The predicted scale-violating NLO Q 2 -dependence is similar to the LO one; 
for x > 0.01 this is also the case for the two rather different input scenarios (3.1), (3.2) 
and (3.1'), (3.2'). In the (x, Q 2 ) region of present data [4-9], A±(x, Q 2 ) increases with Q 2 
for x > 0.01. Therefore, since most present data in the small-x region correspond to small 
values of Q 2 > 0.5 GeV 2 , the determination of g\ (x, Q 2 ) at a larger fixed Q 2 (5 or 10 GeV 2 , 
say) by assuming A^(x, Q 2 ) to be independent of Q 2 , as is commonly done [4-9] (except for 
the recent last reference of 0), is misleading and might lead to an imderestimate of g\ by 
as much as about 20%, in particular in the small-x region x >0.02. (For x < 0.01 the effect 
will be opposite.) The situation is opposite, although less pronounced, for — A™(x, Q 2 ) 
shown in Fig. 3. This implies that \g™ (x, Q 2 ) \ might be overestimated at larger fixed Q 2 by 
assuming v4"(x, Q 2 ), as measured at small Q 2 , to be independent of Q 2 . It is obvious that 
the assumption of approximate scaling for A±(x, Q 2 ) is therefore unwarranted and, in any 
case, theoretically not justified as soon as gluon and sea densities become relevant, due to 
the very different polarized and unpolarized splitting functions (anomalous dimensions) 
in the flavor singlet sector. 

In Fig.4 we compare our NLO results for g?(x,Q 2 ) with EMC, SMC and SLAC- 
E142/E143 data as well as with our LO results which are similar to our original LO 
results [|l| . The reason why the LO results are partly larger by more than about 10% than 
the NLO ones is mainly due to the LO approximation where R = in (2.1). Although 
the agreement between the NLO results and experiment has been significantly improved, 
the EMC If and E143 f§ 'data' at fixed values of Q 2 fall still below our NLO predictions 
in the small-x region. This is partly due to the fact that the original small-x A^-data at 
small Q 2 have been extrapolated [|, |6| to a larger fixed value of Q 2 by assuming A{(x, Q 2 ) 
to be independent of Q 2 . According to the increase of A\ with Q 2 in Fig.3, such an 
assumption underestimates g\ in the small-x region at larger Q 2 . On the contrary, our 
results for gf' do not show such a disagreement in the small- a; region when compared with 
the SMC data |5|, 0] in Figs. 4a and 4b where each data point corresponds to a different 
value of Q 2 since no attempt has been made to extrapolate gi(x, Q 2 ) to a fixed Q 2 from 



12 



the originally measured A±(x, Q 2 ). Our predictions for the NLO parton distributions at 
the input scale Q 2 = ij*%lo m eq.(3.6) with the fit parameters given in (3.7) and (3.8) 
are shown in Fig. 5a; the corresponding LO inputs at Q 2 = [i 2 LO in eq.(3.6) with the fit 
parameters given in (3.7') and (3.8') are shown in Fig.5b. The polarized input densities 
in Figs. 5a and 5b are compared with our reference unpolarized NLO and LO dynamical 
input densities of ref. 0] which satisfy of course the positivity requirement (1.3) as is 
obvious from eq.(3.6). The distributions at Q 2 = 4GeV 2 , as obtained from these inputs 
at Q 2 = yU 2 for the two scenarios considered, are shown in Fig.6. Since not even the 
polarized NLO gluon density 5g(x, Q 2 ) is strongly constrained by present experiments, 
we compare our gluons at Q 2 = 4 GeV 2 in Fig. 7 with the ones which originate from 
imposing extreme inputs at Ql = fi% LO , such as 5g = g, 5g = xg and 5g = 0, instead 
of the one in (3.6) for the 'valence' scenario. The results are similar if these extreme 
gluon-inputs are taken for the 'standard' scenario in (3.6), and the variation of 5g(x,Q 2 ) 
allowed by present experiments is indeed sizeable. This implies, in particular, that the 
(Revolution of gi(x,Q 2 ) below the experimentally accessible x-range is not predictable 
for the time being. 

Finally let us turn to the first moments (total polarizations) Af(Q 2 ) of our polarized 
parton distributions, as defined in (1.1), and the resulting rf' n (Q 2 ) in (3.3). It should be 
recalled that, in contrast to the LO, the first moments of the NLO (anti)quark densities 
do renormalize, i.e. are Q 2 -dependent, due to the non-vanishing of the 2-loop Sj^ 1 in 
(A. 9) and S^g(r] = +1) in (2.10). Let us discuss the two scenarios in turn: 

'standard' scenario: From the input distributions (3.6) together with (3.7), being con- 
strained by (3.1) and (3.2), one infers in NLO 



Q 2 (GeV 2 ) 


Ally 


Ad v 


Aq 


As = As 


Ag 


AS 


r? 


1 1 


^NLO 


0.9181 


-0.3392 


-0.0660 


-0.0660 


0.507 


0.183 


0.1136 


-0.0550 


1 


0.915 


-0.338 


-0.067 


-0.068 


0.961 


0.173 


0.124 


-0.061 


4 


0.914 


-0.338 


-0.068 


-0.068 


1.443 


0.168 


0.128 


-0.064 


10 


0.914 


-0.338 


-0.068 


-0.069 


1.737 


0.166 


0.130 


-0.065 
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Table 1. First moments Af of polarized NLO parton densities Sf(x,Q 2 ) and of 
gi' n (x, Q 2 ) as predicted in the 'standard' scenario. Note that the marginal differ- 
ences for Aq and As indicate the typical amount of dynamical SU(3) / breaking as 
mentioned in footnote 4. 

In LO the input distributions (3.6) together with (3.7') imply 

Au v = 0.9181 , Ad v = -0.3392 , , Aq = As = As = -0.0587 , 
Ag(fi 2 LO ) = 0.362 , Ag(A GeV 2 ) = 1.273 , Ag(10 GeV 2 ) = 1.570 (3.9) 

which result in AS = 0.227. This gives, using eq.(3.5) without the factor a s /ir 
together with (3.1) and (3.2), 

rf = 0.1461 , = -0.0635 (3.10) 

which is similar to our previous LO result |l|. Both our NLO results in Table 1 
and the LO ones in (3.10) are in satisfactory agreement with recent SMC and EMC 
measurements f|, [5|, [7] 

r?(10GeV 2 ) = 0.142±0.008±0.011, r™(10GeV 2 ) = -0.063±0.024±0.013 (3.11) 

as well as with the most recent E143 data @ implying r™(2GeV 2 ) = -0.037 ± 
0.008 ±0.011. 

'valence' scenario: From the input distributions (3.6) together with (3.8), being con- 
strained by (3.1') and (3.2'), one infers in NLO 



Q 2 (GeV 2 ) 
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1 1 
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-0.3392 
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0.496 


0.268 


0.1142 


-0.0544 
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0.915 


-0.338 
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-2.5 x 10~ 3 


0.982 


0.252 
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-0.061 
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0.914 


-0.338 


-0.081 


-3.5 x 10~ 3 


1.494 


0.245 


0.128 


-0.064 


10 


0.914 


-0.338 


-0.081 


-3.8 x 10- 3 


1.807 


0.244 


0.130 


-0.065 



Table 1'. First moments Af of polarized NLO parton densities Sf(x,Q 2 ) and of 
gi' n (x, Q 2 ) as predicted in the 'valence' scenario. 
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In LO the input distributions (3.6) together with (3.8') imply 

Au v = 0.9181 , Ad v = -0.3392 , ,Aq = -0.0712 , As = As = , 

Ag(fi 2 LO ) = 0.372 , Ag(A GeV 2 ) = 1.361 , Ag(10 GeV 2 ) = 1.684 (3.9') 

which result in AE = 0.294. This gives, using eq.(3.4) without the factor a s /ir 
together with (3.1') and (3.2'), 

T? = 0.1456 , = -0.0639 (3.10') 

which is again similar to our previous LO result |IJ. Both our NLO results in Table 
1' and the LO ones in (3.10') compare again well with the experimental results in 
(3.11). 

Apart from the Q 2 -dependent Ag(Q 2 ) in LO and NLO, the Q 2 -dependent first moments of 
NLO (anti)quark densities in Table 1 and 1' should be compared with the Q 2 -independent 
LO results as discussed in the Introduction which, in absolute magnitude, are similar to the 
NLO ones. Although present scarce data obviously cannot uniquely fix the polarized sea 
and gluon densities, our optimal fits favor a sizeable total gluon helicity, Ag(10GeV 2 ) ~ 
1.7, despite the fact that Ag{Q 2 ) decouples from the full first moment Yi{Q 2 ) in (2.13) 
in the MS scheme. 

In both scenarios the Bj0rken sum rule manifestly holds in LO due to our constraints 
(3.1) and (3.1'), and together with the NLO correction due to eq.(3.3) we have 

rW 2 )-rW 2 ) = ^(i-^5) • (3.12) 

It is also interesting to observe that at our low input scales Q 2 = n 2 LO NLO = 0.23, 0.34 GeV 2 
the nucleon's spin is dominantly carried just by the total helicities of quarks and gluons, 
|AS(/i 2 VLO ) + Ag(fi 2 NLO ) ~ 0.6 according to Tables 1 and 1', and |AS + Ag(^i 2 LO ) ~ 0.5 
according to eqs.(3.9) and (3.9'), which implies for the helicity sum rule 

l -= l -A*L{Q 2 ) + Ag{Q 2 ) + L z {Q 2 ) (3.13) 

that L z (ii 2 LO NLO ) — 0. The approximate vanishing of the latter non-perturbative angular 
momentum, being build up from the intrinsic carried by partons, is intuitively expected 
for low (bound-state-like) scales but not for Q 2 3> nlo . 
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4 Summary 



Based on a recent complete NLO calculation |TT| of all spin-dependent two-loop split- 
ting functions 5P-j\x), i,j = q,g (or, equivalently, anomalous dimensions 8^ ) in the 
conventional MS factorization scheme, we have first presented a consistent NLO formu- 
lation of the (Revolution of polarized parton distributions. For calculational purposes 
we have concentrated on (Mellin) n-moments of structure functions where the solutions 
of the NLO evolution equations can be obtained analytically for the parton densities. 
Using these formal results we have performed a quantitative NLO analysis of the longi- 
tudinal spin asymmetry A^ n {x, Q 2 ) and of gi' n (x, Q 2 ), and we have updated our previous 
LO results JTJ. Within the whole relevant x- and Q 2 -region (x > 1CT 3 , Q 2 > 1 GeV 2 ) 
we found a remarkable perturbative stability between LO and NLO results. The scale 
violating Q 2 -dependence of Ai' n (x,Q 2 ) turned out to be similar to the one obtained in 
LO and is non-negligible for (x, Q 2 ) values relevant for present data. The assumption 
of approximate scaling for A±(x, Q 2 ) is therefore unwarranted and theoretically not jus- 
tified. We presented two plausible sets of polarized LO and NLO (MS) parton densities 
Sf(x, Q 2 ) which describe all presently available data very well. In contrast to polarized 
quark and antiquark densities, the gluon density 5g(x, Q 2 ) is rather weakly constrained by 
present data. Our optimal fits, however, favor a rather sizeable total gluon helicity, e.g., 
Ag(Q 2 = 10 GeV 2 ) ~ 1.7. It should be reemphasized that only processes where 5g occurs 
directly already in LO (with no 5q and Sq contributions present) appear to be the most 
promising sources for measuring Sg(x, Q 2 ). This is the case for 7* (7)^0 — > cc responsible 
for open charm or J/\l/ production (see, e.g., ref. [15|]). Our results demonstrate the com- 
patibility of our restrictive radiative model, cf. eq.(1.3), down to Q 2 = ^nlo = 

0.34 GeV 

and to Q 2 = n 2 LO = 0.23 GeV 2 , with present measurements of deep inelastic spin asym- 
metries and structure functions. 

A Fortran package containing our optimally fitted 'standard' and 'valence' NLO 
(MS) as well as LO distributions can be obtained by electronic mail from 
stratmann@het.physik.uni-dortmund.de or vogelsang@v2.rl.ac.uk. 
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Appendix 



The spin-dependent LO anomalous dimensions (splitting functions 1 ) have been originally 
calculated in |34|, |35| and are given by 



Iqq 



4C F 
-87) 

4^4 



2S 1 (n) 
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n(n + 1) 
2S 1 {n) - 



-4CV 
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11 

IF 



(A.l) 



n(n + 1) 

where Cp = 4/3, Ca = 3 and 7) = //2, with / being the number of active flavors (/ = 3 
has been used when calculating 5jij). Note that dj^s = 8 n ff^ n = 7^ n where the latter 
quantity refers to the spin-averaged (unpolarized) anomalous dimension. Furthermore, 
for the first n = 1 moment we have S^yjj^ 1 = 5^yJ l = as a consequence of helicity 
conservation at the quark-gluon vertex. 

The spin-dependent NLO (MS) two-loop flavor non-singlet anomalous dimensions 
£7^5 (77), required in (2.10) for the evolution of 5q^ s V= ±(Q 2 ), are the same as found for 
the spin-averaged case, S^^siv) = I^Nsiv) with 7^ n (^ = ±1) being given by eq.(B.18) 
of |||. Note that S^iv = +1) governs the evolution of the NS combinations 5q — 5q, 
while S^girj = —1) refers to the combinations 5q + 5q appearing in the NS expressions 
(2.11). The NLO flavor singlet anomalous dimensions SjQ in the MS scheme are as 



follows 01: 



Iqq 
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with 7tvs' 1 (?7 = — 1) being again given by eq.(B.18) of J36| and 5 

n 4 + 2n 3 + 2n 2 + 5n + 2 



h [ pl n m = iec F T f 



iqg 



«J 7 (i) 

'99 



where 



8CpTf 



n 



n 3 (n + l) 3 
L-(£ a (n)-5?(n))+4 n ~ 1 



(A.3) 
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1. This latter property has, howe 
positivity, eq.(1.3), of our resulting NLO parton distributions. 
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with G(n) = ^ - i> (f) and 77 = (-) n -> ±1 for = ±1) and 77 — > — 1 for 

the flavor singlet anomalous dimensions (evolutions). The analytic continuations in n, 
required for the Mellin inversion of these sums to Bj0rken-x space, are well known [ TB\ . 



It should be noted that the original results for 5^ have been presented [JTTJ] in 
terms of multiple sums, denoted by Sk{n), Sk t i(n) and Skj(n), which cannot be directly 
analytically continued in n. The following relations have been used in order to arrive at 
(A.4)-(A.6): 
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where for the latter sum we have used the identity 
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in order to relate §i >2 to the expressions in (A.7). 



Finally, the first n = 1 moments of the q — > q(q) and g — > g(g) transitions reduce to 
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Figure Captions 



Fig.l Comparison of our NLO and LO results for A±(x, Q 2 ) as obtained from the fitted 
inputs at Q 2 = fi% LO lo f° r the 'standard' (eqs.(3.7) and (3.7')) scenario with present 
data [4-9]. The Q 2 values adopted here correspond to the different values quoted in 
[4-9] for each data point starting at Q 2 ^ 1 GeV 2 at the lowest available x-bin. The 
results in the 'valence' scenario are indistinguishable from the ones shown. 

Fig.2 The Q 2 -dependence of A p { d (x,Q 2 ), as predicted by the NLO and LO QCD evolu- 
tions at various fixed values of x, compared with recent SLAC-E143 || (solid circles) 
and SMC data |5|, ^] (open circles). 

Fig. 3 The Q 2 -dependence of A\' n {x, Q 2 ) as predicted by the NLO and LO QCD evolu- 
tions at various fixed values of x. 

Fig. 4a Comparison of our 'standard' and 'valence' NLO and LO results with the data [4- 
9] for g p ' d (x, Q 2 ). The SMC data correspond to different Q 2 > 1 GeV 2 for x > 0.005, 
as do the theoretical results. 

Fig. 4b Same as in Fig. 4a but for g™(x,Q 2 ). The E142 and E143 data || §] correspond 
to an average (Q 2 ) = 2 and 3 GeV 2 , respectively, and the theoretical predictions 
correspond to a fixed Q 2 = 3 GeV 2 . 

Fig. 5a Comparison of our fitted 'standard' and 'valence' input NLO (MS) densities in 
eqs.(3.7) and (3.8) with the unpolarized dynamical input densities of ref. [14]. 



Fig. 5b The same as in Fig. 5a but for the LO input densities according to eqs.(3.7') and 
(3.8'). 

Fig. 6 The polarized LO and NLO (MS) densities at Q 2 = 4 GeV 2 , as obtained from the 
input densities at Q 2 = nfjLo lo * n Figs.5a,b. In the 'standard' scenario, 5s coincides 
with the curves shown for Sq in NLO and LO due to the SU(3) f symmetric input 
which is only marginally broken in NLO for Q 2 > [i 2 NLO - 

Fig. 7 The experimentally allowed range of NLO polarized gluon densities at Q 2 = 4 GeV 2 
for the 'valence' scenario with differently chosen 8g(x, n 2 NLO ) inputs. The 'fitted 5g' 
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curve is identical to the one in Fig.6 and corresponds to 5g(x, fj, 2 NLO ) in eq.(3.8). 
Very similar results are obtained if 5g(x, fi 2 NLO ) is varied accordingly within the 
'standard' scenario as well as in a LO analysis (see, e.g., ref.|Q). 
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